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ABSTRACT 

In this paper we present a deep and homogeneous i-band selected muhi-waveband 
catalogue in the COSMOS field covering an area of about 0.7n°. Our catalogue with 
a formal 50% completeness limit for point sources of i ~ 26.7 comprises about 290 000 
galaxies with information in 8 passbands. We combine publicly available u, B, V, 
r, i, z, and K data with proprietary imaging in H band. We discuss in detail the 
observations, the data reduction, and the photometric properties of the H-band data. 
We estimate photometric redshifts for all the galaxies in the catalogue. A comparison 
with 162 spectroscopic redshifts in the redshift range < z < 3 shows that the 
achieved accuracy of the photometric redshifts is Az/(zspec + l) < 0.035 with only 
~ 2% outliers. We derive absolute UV magnitudes and investigate the evolution of 
the luminosity function evaluated in the restframe UV (1500 A). There is a good 
agreement between the LFs derived here and the LFs derived in the FORS Deep 
Field. We see a similar brightening of M* and a decrease of 4>* with redshift. The 
catalogue including the photometric redshift information is made publicly available. 
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1 INTRODUCTION 

In the last decade our knowledge about the evolution of 
global galaxy properties over a large redshift range has 
improved considerably. The 2dF Galaxy Redshift Survey 
(2dFGRS; ,Collcss et aLll200ll). th e Sloan Digital Sky Sur- 
yev (SDSS:lstoughton et al.ll2002l ). and the 2MASS survey 
iJarrett et al.ll200(]| ) have provided very large local galaxy 
samples with spectroscopic and/or photometric information 
in various passbands. Thanks to these data sets we are now 
able to assess very accurate local [z ~ 0.1) reference points 
for many galaxy evolution measurements like the luminosity 
function, the star formation activity, the spatial clustering 
of galaxies, the stellar population, the morphology, etc. 

In the redshift range between 0.2 < z < 1 pioneer- 
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f Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA), operated jointly by the Max-Planck- 
Institut fiir Astronomic, Heidelberg, and the Instituto de As- 
trofisica de Andalucia (CSIC). 



ing work has been do ne in the contex t of the Canada 
Fran ce Redshift Surv ey (|Lilly et al.lll995l ). the Autofib sur- 
vey l|Ellis et al.ll 19961 ) and in th e Canadian Network for Ob- 
servational Cosmology survey (|Yee et al.lll996 l. They pro- 
vide accurate distances and absolute luminosities by spec- 
troscopic foUowup of optically selected galaxies, thus being 
able to probe ba sic properties of gal axy evolution. Moreover 
the K20 -survev (ICim atti ct al. 200 ^ as well as t he MUNICS 
survey (|Drorv et al., ,2001a : Fculner et al.l bOOd ) extend the 
analysis into the near infrared regime (for 0.2 < z < 1.5). 

An important step towards probing the galaxy proper- 
ties also in the high redshift regime aroun d z ~ 3 and z ^ 4 
was th e work of lSteidel Sz HamiltonI l|l993l ) and lSteidel et al] 
(119961 ). They used colour selection to discriminate between 
low and high redshift galaxies (see lGiavaliscoll2002l . for a re- 
view). The so-called Lyman-break galaxies (LBGs, mainly 
starburst galaxies at high redshift) are selected by means 
of important features in the UV spectrum of star-forming 
galaxies. 

The next milestones in pushing the limiting magni- 
tude for detectable galaxies to fainter and fainter limits 
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were the space based Hubble Deep Field North (HDFN; 
IWilliams et"aLl Il996al) and Hubble Deep Field S outh 
fHDFS: I Williams et all l2000l : ICasertano etaP l200d ) (see 
[Ferguson et al. 2000l . for a review). Although of a limited 
field of view of about 5D' only, the depth of the HDFs al- 
lowed the detection of galaxies up to a redshift of 5 and even 
beyond. 

In the past years the space based HDFs were supple- 
mented by many more multi-band p hotometric surveys lik e 
the NTT SUSI deep Field (NDF; lArnouts et all Il999all 
the C handra Deep Field South (CDFS: lArnouts et al.l 
200181). the William Herschel Deep Fie ld (WHDF; 



McCracken et al ] |2000al : [Metcalfe et al.l 



Deep Field/Survey fSDF: lMaihara et al 
I2OO3I). the COMBO -17 survev ijWolfet al 



2001al:l0uchi et all 
""^OOSl), FIRES 



jLabbe et all l2003t ). the FORS Deep Field (FDF, 
iHeidt et aLl bOOSl ) the Great Observator ies Origins Deep 
Survey (GOODS: iGiavalisco et all l2004bl) the Ultra Deep 
Field ( UDF and UDF- Parallel ACS fields: I Giavalisco et al.) 
l2004al : iBunker et aLll2004l: iBouwens et al.ll2004l ) the VIR- 
MOS deep survev (|Le Fevre et al.ll2004l')rGEM S (Rix et al' 
I2OO4I '). the Keck Deep Fields (|Sawicki fc Thom pson 2005), 
and the Multiwave l ength Survey by Y ale-Chile (MUSYC; 
iGawiser et al]|2006l : lOuadri et al.ll2007l ). 

With the advent of all these deep multi-band pho- 
tometric surveys the photometric redshift technique (es- 
sentially a generalisation of the drop-out technique) can 
be used to identify high-redshift galaxies. Photomet- 
ric redshifts are often determined by means of tem- 
plate matching algorithm that applies Bayesian statis- 
tics and uses semi-empirical temp late spectra matched to 
broad-band photom e try (see also Baum 19621: Kool 1985|: 
Brunner et al. I 1 19991: li^nandez-Soto et all Il999r iBemtea 



2000| : iBender et all bOOll : iLe Borgne fc Rocca-Volmerangd 
20021 : iFirth et allbOOd ). Redshifts of galaxies that are sev- 



eral magnitudes fainter than typical spectroscopic limits can 
be determined reliably with an accuracy of Az/(zspec + 1) 
of 0.02 to 0.1. 

In this context the COSMOS survey ("Scoville et alj 

l|2007h : see also http://www.astro.caltech.edu/cosmos/ for 
an overview) combines deep to very-deep multi-waveband 
information in order to extend the analysis of deep pencil 
beam surveys to a much bigger volume, thus being able to 
drastically increase the statistics and detect also very rare 
objects. For this, the survey covers an area of about 2D° 
with imaging by space-based telescopes (Hubble, Spitzer, 
GALEX, XMM, Chandra) as well as large ground based tele- 
scopes (Subaru, VLA, ESO-VLT, UKIRT, NOAO, CFHT, 
and others). 

In this paper we combine publicly available u, B, V, 
r, i, z, and K COSMOS data with proprietary imaging in 
the H band to derive a homogeneous multi-waveband cat- 
alogue suitable for deriving accurate photometric redshifts. 
In Section [5] we give an overview of the near-infrared (NIR) 
data acquisition and we describe our 2-pass data reduc- 
tion pipeline used to derive optimally (in terms of signal- 
to-noise for faint sources) stacked images in Section [3] We 
also present NIR galaxy number counts and compare them 
with the literature. 

In Section |3] we present the deep multi-waveband i-band 
selected catalogue and discuss its properties, whereas the 
data reduction of the spectroscopic redshifts is described in 



Section [5l In Section [6] we present the photometric redshift 
catalogue, discuss the accuracy of the latter and show the 
redshift distribution of the galaxies. In Section [7] we derive 
the redshift evolution of the restframe UV luminosity func- 
tion and luminosity density at 1500 A from our i-selected 
catalogue before we summarise our findings in Section IS] 

We use AB magnitudes and adopt a A cosmol- 
ogy throughout the paper with Q.m = 0.3, Q,a ~ 0.7, and 
Ho = 70kms-'Mpc-^ 



2001alV the Subaru 2 NIR OBSERVATIONS 



2.1 Field layout 

Our observing strategy was designed to follow-up the pub- 
licly available COSMOS observations with proprietary imag- 
ing in H band. The whole area is covered by 25 point- 
ings (15.4'xl5.4' each) to a depth adequate to the public 
NIR data. The field layout and nomenclature is shown in 
FigEl The programme was carried out as a joint effort be- 
tween three large extragalactic survey projects currently be- 
ing pursued at the Centro Astronomic o Hispano Alema n 
(CAHA), on C alar Alto: ALHAMBRA (iMo les et al.' 1 20051 ') . 
MANOS-wide (|Roser et al.ll2004 IZatloukal et al. 200'^) and 
MUNICS-Deep (Goranova et al., in prep.). In addition, as a 
part of the MUNICS-Deep project, in one of the pointings 
(01c) we have already collected deep Js- and K'-band data. 



2.2 Data acquisition 

The observations in H, Js and K' bands were carried out 
using the NIR wide- field imager OMEGA200CE], operating at 
the prime focus of the CAHA 3.5m telescope. OMEGA2000 
is equipped with a HAWAII-2 HgCdTe 2048x2048 array. The 
instrument pixel scale is 0.45"per pixel, providing a field-of- 
view of 15.4'x 15.4'. 

Here we present the H-band observations in 15 point- 
ings (see FigUI, collected during 11 nights spanned over 3 
observing campaigns: December 2004, February 2005, and 
March/April 2005. Each pointing was observed for a total 
of at least 3 ksec. The individual exposure times were 3 sec, 
on-chip co-added to produce single frames of 60 sec each. 
Except for 05b, at least 50 such frames (depending on the 
weather conditions) were observed per pointing. 

In addition, observations in Js and K' bands in 01c were 
collected during 7 night in the observing campaigns Novem- 
ber 2003 and February 2006. The total exposure time in Js- 
band data was 8.2 ksec with the same individual exposure 
time scheme as for the H band. The K'-band observations 
have a total of 7.7 ksec exposure time with individual expo- 
sures of 2 sec, co-added internally into single frames of 30 
sec each. 

All observations were done using dithering pattern con- 
sisting of typically 20 positions shifted with respect to one 
another by 20". The consecutive dithering sequences were 
repeated with the same pattern but with an offset in the 
origin. 



1 http://w3.caha.es/CAHA/Instruments/O2000/index2.html 
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Table 1. The 15 COSMOS patches. The table gives the filter, the 
field name, the coordinates (right ascension a and declination S) 
for the equinox 2000, the number of frames, and the total exposure 
time. 

filter field a 5 frames exp. time 

(2000) (2000) [ksec] 



d) Rejecting cosmic rays, 

e) Aligning the images, and 

f) Stacking the images 

• Second pass: 

g) Constructing the object mask from the stacked image 

h) Masking all objects in the raw frames using the 
object-mask, repeating pass 1 (step a and b) with masked 
objects and proceed as in the 1^^ pass for step c, d, e, and 
f until final stacking. 

The reasons for using this 2-pass reduction pipeline are 
the following: Due to severe problems with our dome and 
twilight sky flats (see Goranova et al. in prep.), we derive the 
flat-fields from the science frames themselves. In such a case, 
however, it is mandatory to to ensure that the final flat-field 
is free of object residuals. Due to the very high noise level in 
the raw NIR data, this is possible only if all objects (both the 
bright and the faint ones) are masked out properly before 
deriving the flat-field, which we could not achieve with any 
other objects rejection algorithms {k — a clipping, min-max 
rejection, etc.) we have applied. An example is presented in 
Fig. [1] where we show the ratio between the flrst pass (aver- 
aging after min-max rejection) and second pass (averaging 
after min-max rejection and object masking) flat-fields. The 
black regions clearly show the object residuals in the flat- 
field used for the first pass. Although the relative systematic 
error is in the order of only 0.1%, it makes a significant ef- 
fect in the final stacked image. The reason for this is the 
very high sky level of the NIR images. In our case the sky 
level A'^ of a stacked image is at least in the order of 9 mil- 
lion photons. This implies, that the Poissonian error (\/iV) 
is about 3000 photons. If we now introduce a systematic er- 
ror of 0.1%, we end up with an error of 9000 photons. It is 
therefore clear that we introduce a systematic error of about 
3(7. For faint objects, which are detected on the level of a few 
(J of the sky background, the systematic offset is in the same 
order as their total flux. In other words, if we introduce a 
systematic flat-field error of 0.1%, the introduced systematic 
error is on the same level as the statistical Poissonian error 
for a sky level of 1 million photons. 

Another problem of NIR data reduction is the short- 
scale time variability of the sky. We subtract the sky using 
the so-called 'running' sky technique, i.e we determine the 
sky for each frame individually. For this we take a series of 
usually 5 frames grouped in time before and after (but ex- 
cluding) the frame we are subtracting the sky from and we 
average those using min-max rejection and object masking. 
This works very well as the sky-illumination does not change 
much within this short time. 

The variability of the sky on the other hand may change 
the overall shape of the flat-field (derived before sky sub- 
traction). In order to avoid correcting this additive illumi- 
nation by the multiplicative flat-field, we use the following 
approach. First we derive flat-fields for each night by av- 
eraging all frames from that night using min-max rejection 
and object masking. Then we average the normalised flat- 
fields from all nights in order to derive a master flat-field. 
This master flat-field should be the best approximation to 
the real shape of the overall sky-flat. Then we fit a 4*''-order 
Chebyshev polynomial to the master flat-field in order to get 
the shape of the master flat-field and eliminate the pixel-to- 



01c 


10:01:26.00 


-1-02:26:41.0 


93 


5.58 


Old 


09:59:31.24 


-f02:26:41.0 


81 


4.86 


02c 


10:00:28.61 


-1-02:26:41.0 


84 


5.04 


02d 


09:58:33.85 


-1-02:26:41.0 


100 


6.00 


03b 


09:59:31.23 


-1-02:41:01.0 


50 


3.00 


03c 


10:01:25.99 


-1-02:12:21.1 


50 


3.00 


03d 


09:59:31.25 


-f02:12:21.1 


50 


3.00 


04b 


09:58:33.84 


-1-02:41:00.9 


74 


4.44 


04c 


10:00:28.60 


-1-02:12:21.0 


49 


2.94 


04d 


09:58:33.86 


-f02:12:21.0 


50 


3.00 


05a 


10:01:25.98 


-f01:58:01.1 


49 


2.94 


05b 


09:59:31.26 


-1-01:58:01.1 


24 


1.44 


06a 


10:00:28.62 


4-01:58:01.0 


70 


4.20 


07a 


10:01:25.98 


-1-01:43:41.1 


61 


3.66 


08a 


10:00:28.63 


-1-01:43:41.1 


48 


2.88 


01c 


10:01:26.00 


-1-02:26:41.0 


136 


8.16 


01c 


10:01:26.00 


-1-02:26:41.0 


258 


7.74 



The observing log presenting filters, pointings nomen- 
clature, coordinates, number of frames, and the total expo- 
sure time per pointing are listed in Table [T] . 



3 NIR IMAGING DATA REDUCTION 

In this section we present the NIR data reduction in the Js, 
H, and K' bands. We develop and describe a 2-pass reduc- 
tion pipeline optimised for reducing and stacking images of 
different quality to get an optimal signal-to-noise ratio for 
faint (sky dominated) pointlike objects. We extract Js, H, 
and K' selected catalogues and show the accuracy of our 
photometric calibration. Based on these catalogues we com- 
pare the number counts in the three NIR filters with number 
counts taken from the literature. 

3.1 Basic reduction 

For the basic reduction we use our own reduction pipeline 
(see Goranova et al., in prep, for a detailed description) 
based on the IRAF0 external package XDIMSU^fl. We use a 
2-pass reduction implementing object-masking for both flat- 
field and sky determination. The two passes are as follows; 

• First pass: 

a) Constructing the 1^^ science-flat, 

b) Subtracting the sky, 

c) Masking bad pixels and bad regions, 

^ "IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation." 

^ Experimental Deep Infrared Mosaicing Software; XDIMSUM is 
a variant of the DIMSUM package developed by P. Eisenhardt, M. 
Dickins on, S.A. Stanford, and J. Ward. F. Valdes (IRAF group); 
see also [Stanford et al. I lll995l '). 
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Figure 3. H-band magnitude differences (black dots) between our COSMOS catalogue and 2MASS in 2 fields (03b and Olc). The red 
symbols show the average in the different magnitude bins. Please note that for these plots we used Vega magnitudes and not AB, in 
order to avoid adding any systematics by transforming the 2MASS magnitudes to AB system. The error in the zeropoint determination 
can be found in Table [2] 




Figure 1. Ratio of first and second pass flat-field: The first pass 
flat-field has been divided by the second pass flat-field. Dark re- 
gions (in the order of 0.1%) represent object residuals in the fiat- 
fleld after the first pass. See text for details. 

pixel variations. Finally, we bend every single flat-field to 
the shape of the Chebyshev flat-field. The corrections ap- 
plied to the single flat-fields are in the order of 1-2% (peak 
to peak). 



Table 2. Field characteristics. The first column gives the filter 
band, the second the field nomenclature as used in Fig. |4] In the 
third column, the number of detected object is reported (see text 
for details). The 50% completeness limit follows in the fourth 
column, the FWHM seeing value in the fifth, and the error of the 
photometric zero point in the last column. 



field 


objects 


50% CL. 


seeing 


zeropoint 






[mag] 


["] 


[mag] 


Olc 


2744 


21.76 


1.13 


0.006 


Old 


4013 


22.38 


0.82 


0.005 


02c 


3483 


22.24 


0.92 


0.007 


02d 


3364 


22.40 


0.98 


0.008 


03b 


2996 


21.82 


1.14 


0.006 


03c 


2020 


21.27 


1.41 


0.007 


03d 


2765 


21.87 


1.24 


0.007 


04b 


2846 


22.25 


1.00 


0.006 


04c 


2567 


21.68 


1.09 


0.006 


04d 


2348 


21.52 


1.39 


0.008 


05a 


1815 


21.16 


1.07 


0.008 


05b 


1102 


20.60 


1.56 


0.006 


06a 


2981 


22.16 


0.83 


0.008 


07a 


1940 


21.52 


0.80 


0.007 


08a 


2618 


21.83 


0.94 


0.019 


Olc 


4005 


22.67 


0.91 


0.009 


Olc 


3174 


21.76 


1.05 


0.009 



Since our final goal is the combination of our H-band 
data with public data to construct a multi-wavelength cat- 
alogue, we want to astrometrically align our H-band im- 
ages to the publicly available images in the COSMOS field 
(jScoviUe et all l2007l ). Moreover, as we intend to use the 
same procedure as in the FDF in constructing the catalogue 
l|Heidt et al.ll2003l ). i.e. using SExtractor in the double im- 
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Figure 4. Field nomenclature of the 15 NIR patches together with an excerpt of the H-band detected catalogues (green, black and white 
for clarity) on top of the HST COSMOS image. Each catalogue covers an area of 205. 44n'. 



age mode, the data sets in the diflFerent passbands must also 
be aligned by pixel. On the other hand, we do not want to 
shift the NIR images twice. This would smooth the images 
and destroy the signal of very faint objects. As it is also 
very challenging to shift single NIR images on top of deep 
optical images (where many bright starlike objects are e.g. 
saturated), we use the following approach: 
First we compute the astrometric transformation between 
our stacked NIR image and the reference optical frame (Sub- 
aru z-band). This allows us to rely on hundreds of objects 
when computing the astrometric solution instead of only a 



few tenths of objects visible in a single NIR image. Knowing 
the relative shifts of the single images to the stacked XDIM- 
SUM image and the full astrometric solution of the latter 
(with respect to the Subaru z-band), we are able to calcu- 
late the astrometric solution also for each single, unshifted 
pre-reduced image. Restarting from these images, we do the 
alignment, the distortion correction, the regridding and the 
stacking of the individual images (with optimal signal-to- 
noise ratio, see below) in one step. 
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15 80 35 15 80 35 



Figure 6. Left panel: Js-band galaxy number counts of patch 01c (red dots, not corrected for incompleteness) as compared to the 
literature (blue dots ) . The vertical dotted line ind i cates the 50% comp leteness limit. The literature number counts are taken from 
ISaracco etal] lll999l . 1200 ll) : iMaihara et al.l l l200ld) : iTeplitz etal] l ll999l ). Right panel: K'-band galaxy number counts of patch 01c 
(red dots, not corrected for incompleteness) as compared to the literature (blue dots ). The vertical dott e d line indicates the 50% 
completeness limi t. The liter ature number counts are taken from Gardner ct al. (1993); Soifcr ct al.' ('1994'); Diorgovsk^_et_al1 
Glazebrook et al.l([T9 95l: McLe od et"all lll995l'):lGa rdner et a l. (1996); Huang ct al. (199 7); Moustakas ct al. (1997'): Saracco et al.l 



199£); 



1997); 



Bershadv et al.l(ll998 ); Szokoly ct al. (1998); McCrackcn ct al. (2000b); Vais anen et al.l (20001: Drory ct al. ( 2001bl)- lH'uane et al (200M: 
Kummel fc WagneJ ( 12001.) : iMaihara et al, (, 2001b ): ,Martinii t2001h ): iSaracco et al.l ( l200ll' ): .Crist6bal-Hornillos et all (l2003 l ) : lMinowa et all 



(|2005l V The error bars show the la Poissonian errors. 



3.2 Image stacking with optimal S/N for point 
sources 

As the single exposures were not taken under the same ob- 
serving conditions, the sky levels, the seeing, and the zero- 
points can substantially differ from frame to frame. On the 
other hand this gives one the possibility to stack the sin- 
gle images with weighting factors in order to achieve a final 
combined frame with optimal signal-to- noise ratio (S /N) for 
point sources. 

We calculated a weight a to be applied to an individual im- 
age following the general Ansatz for two images (denoted by 
index 1 and index 2): 



Stot 



Si + 012S2 



Ntot = y^iV2 + (a2iV2)2 



(1) 



(2) 



where Stot and Ntot is the signal and the noise of the com- 
bined image. This transforms to: 

Stot _ fi + Q2/2 

Ntot ^ V(/i + hiafn) + ai{f2 + hia'in) 

where /i and /2 are the fluxes of an object without sky (sig- 
nal), hi and /12 are the sky values and ai and (T2 correspond 
to the seeing in the two frames. 02 is the weighting factor to 
be applied to frame 2. It is than straight forward to compute 
the value of Q2 for which Stot /Ntot is maximised: 

" Nto 



da2 



= 



ce2 — 



/2(/i + feiaiV) 

M.f2 + h2al7T) 



(3) 



For bright objects which are not dominated by the sky noise 
(/ 3> h) Equation ([3| approaches the limit 



Q2 = 1, 



(4) 



whereas for faint, sky dominated objects (/ <C h) Equa- 
tion ((3]) transforms into 



Q2 



/2 ■ (hia!) 



{h2al) 



(5) 



As the overwhelming majority of the objects are very faint 
point sources and therefore dominated by the sky noise. 
Equation ((5} is used to derive the weighting-parameter ct. 

To reduce errors when determining the weighting- 
parameter Of, the fluxes /1/2 are derived from a bright not 
saturated star, the sky levels h-ij2 correspond to the mod^B 
of the image and the seeings (T1/2 have been calculated from 
the median seeing of a few stars. 

The weighting factor for the frame (randomly cho- 
sen) was set to unity. The factors for all other (N — 1) im- 
ages were then derived following Equation ((5} relative to this 
image. The flnal stacked frame /sum can then be calculated 
according to: 



(6) 



most frequent value in the pixel histogram of the image 
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Figure 5. Galaxy number counts in the H-band (not cor- 
rected for incompleteness) as compared to the literature (blue 
dots). The black dots represent the number counts of the single 
patches, whereas the red dots show the average number counts 
of all fields. The vertical dotted lines indicates the 50% com- 
pleteness limits of the shallowest (05a) and the dee pest (02d) 
patch es. The literatur e numb e r counts are taken from I Ya ii^t^^l. 
1 1998l)j iTeplitz et aTl h99§): 'Thompson et al.' (1999'); Martini 
l2001a^ -IChen et al.l ^20 021: Mov ct al. (2003); Firth ct al. (200f|; 
iMetcalfe et al. I l l200d) . The error bars show the la Poissonian er- 
rors. 



Figure 2. Different steps of the reduction pipeline: A single raw 
H-band image (upper panel), the same image after pre-reduction 
(middle panel), and a stacked image after the two-pass data re- 
duction and stacking with optimal S/N ratio (lower panel). 



where the index i denotes a single image, ai is the weighting 
factor according to Equation ((5|; h is the sky-subtracted 
single image; Bi is the bad pixel maslfl and fi is the flux 
derived from the bright star (used to derive a;; see also 
Equation [2} . 

Finally, we did the image alignment, the distortion cor- 
rection, the regridding and the stacking of the individual im- 
ages with optimal S/N in one step by using standard IRAF 
routines as well as SWarp (|Bertin. E.II2003I ). The regridding 
has been done to the native pixel scale of the Subaru tele- 
scope (0.2" per pixel). Please note that we do not interpo- 
late any bad pixel or bad region in the single images, but set 
them to zero during the stacking procedure. Missing flux in 
bad regions is taken into account by using Equation (|6]). 
Since a diflerent number of dithered frames contributed to 
each pixel in the co-added images (producing a position- 
dependent noise pattern) a combined weight map for each 
frame was constructed. The latter was used during source 
detection and photometry procedure to properly account for 
the position-dependent noise level. 

In Fig. [2] we illustrate different steps of our reduction 
pipeline: a single raw H-band image, the same image after 
pre-reduction, and a stacked image after the two-pass data 
reduction and stacking with optimal S/N ratio. There are 



zero for a bad pixel, unity otherwise 
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practically no haloes around very bright objects, e.g. due 
to object residuals in the flat-field or sky subtraction prob- 
lems / residuals. 

3.3 Photometric calibration 

The absolute photom etric calibration is based on the 2MASS 
l|jarrett et al.|[200ol ) catalogue. First we cross-correlate the 
objects in all our reduced patches with sources in the 2MASS 
catalogue. We exclude all objects with possible problems in 
either one of the catalogues relying on the quality flags of 
2MASS and SExtractor quality flags (only objects with a 
flag of ^ 3 are considered). An error weighted flt between 
these objects (20 to 45, depending on the patch) then deter- 
mines the zeropoint as well as its error. In Fig. [3] we show 
typical magnitude differences between our COSMOS cata- 
logue and 2MASS. Please note that for these plots we used 
Vega magnitudes and not the AB magnitude systems (we 
did not want to add any systematic by changing the 2MASS 
magnitude system) . The accuracy of the zeropoint in the dif- 
ferent flelds is in the order of 0.01 magnitudes (derived from 
the error weighted flt between the COSMOS and 2MASS 
objects; see above) and can be found in Tabled 

3.4 NIR catalogues 

Based on the stacked images we derived 15 H-band selected 
catalogues (see Fig. |3| as well as 1 J-band and 1 K-band 
selected catalog ue (in fleld 01c, only). F or this purpose we 
run SExtractor (|Bertin fc Arnout3l 19961 ) with the detection 
threshold t = 2 (minimum signal-to-noise ratio of a pixel 
to be regarded as a detection) and n = 3 (number of con- 
tiguous pixels exceeding this threshold). Depending on the 
deptljfl of the different patches, we detect between 2000 and 
4000 objects (excluding patch 05b were we have only half 
of the minimum exposure time we wanted to achieve). The 
false detection rate are in the order of one percent or less 
(detected on the inverted image). Only around extremely 
bright and saturated objects (three in patch 02d and one in 
patch 03c) the false detection rate increases. As the regions 
around these four objects are not taken into account, we get 
an overall false detection rate in the order of 1 percent. 
Because the depth of the patches decreases towards the bor- 
ders, we limited our analysis to the inner fleld. The signal- 
to-noise ratio in this 'deep' region is at least 90% of the 
best S/N in every patch. This prevents a possible bias of 
the photometric redshifts due to a not completely homoge- 
neous data set. The single patches, each covering an area of 
205. 44D' together with the field nomenclature are shown in 
Fig. U In total we detected about 40 000 H-selected objects 
over an area of about 0.85D° as well as about 4000 (3000) J 
(K) selected galaxies over an area of 205. 44D'. 

3.5 NIR number counts 

The galaxy number counts can be used to check the calibra- 
tion of the data set, to detect possible galaxy over or under- 
densities of a field as well as to determine the approximate 

^ as a result of the slightly varying seeing and total exposure 
times, see Table [2] 



Table 4. COSMOS (PDF) field characteristics for a SExtractor 
detection threshold of t = 2.5 (t = 1.7) and n = 3 (n = 3) 
contiguous pixels. 



Filter 


50 % CL. 


50 % CL. 


seeing 




COSMOS 


FDF 


["] 


UCFHT 


25.6 


26.5 (U) 


0.90 


^Subaru 


27.7 


27.6 (B) 


0.95 


^Subaru 


26.5 


26.9 (g) 


1.30 


^Subaru 


26.8 


26.9 (R) 


1.05 


^Subaru 


26.7 


26.8 (1) 


0.95 


'^Subaru 


25.1 


25.8 (z) 


1.15 




21.2 


22.6 (Ks) 


1.28 



depth of the data. We did not put much effort in star-galaxy 
separation at the faint end (as we did for the i-selected cata- 
logue, see Sect.|4|, where the galaxies dominate the counts. 
At the bright end, where SExtractor is able to disentangle a 
stellar and a galaxy profile, we used the star-galaxy classifier 
to eliminate obvious stellar objects. We present in Fig.[5]the 
H-band number counts of all patches and in Fig. [6] the Js 
and K' band number counts of patch 01c. Although the sin- 
gle patches show a scatter in the H-band number counts at 
the very bright end, there is a good to very good agreement 
between the literature number counts and our mean number 
counts ( red dots) up to the 5 % completeness limit for point 
sources (|Snigulaetal .1120021 ). Comparing the Js and K' band 
number counts with data taken from the literature (Fig. O 
also shows a relatively good agreement, although patch 01c 
seems to be overdense with respect to most of the literature 
values. This is also true for the H-band number counts in 
01c if compared to the other patches or to the literature, 
indicating an overdensity in this specific patch. Neverthe- 
less, our results are compatible within the error bars at a la 
level. 

Please note that the 50 % completeness levels for point 
sources (vertical red lines in Fig. [S] and Fig.|Sl) coincide very 
nicely with the faint-end region where the number counts 
start do drop. The values of the NIR galaxy number counts 
can be found in Table O 



4 THE I-SELECTED CATALOGUES 

Based on the publicly available optical and NIR data of the 
COSMOS fiel4!| we build a Subaru i-band detected galaxy 
catalogue in 12 of our 15 patches. Because of the relatively 
bad H-band seeing (> 1.3"), we exclude the patches 03c, 04d 
as well as 05b. We decided to use the ground based i-band 
for source detection mainly for two reasons: First, at the 
time when we started the photometric redshift determina- 
tions the HST coverage of our reduced H-band patches was 
not yet completed. Second, including space based data with 
their superb PSF would need a different technique (with 
HST images in various wavebands) to optimally construct 
the multi-waveband catalogue (see e.g. Grazian ct al. ( 20061) 
for the GOODS-MUSIC sample). 

We use a very similar approach as for the FDF l|Heidt et al.l 

The data were taken from: 
|http://irsa.ipac. caltech.edu/data/COSMOS/ 1 
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Figure 7. Left panel: Galaxy number counts in the i-band (not corrected for incompleteness) as compared to the literature. The black 
dots represent the number counts of the single COSMOS patches, whereas the red dots show the average number counts o f all fie l ds. The 
vertic al dotted line indica t es the 50% co mpleteness limit. The li t eratur e number counts are tak e n from [Hal^fc Mackav ('l984'') :lTvson| 



198g ): iLillv et al.l lll99j): ICasertano et a l. (199 5.): IPriver et al.l lll995l): IWiUiams et al.] lll996d): iHuang et aTT(il 998): Postm an et al, 
1998l i: lArnouts et al.l l ll999bh : iMetcalfe'et al.i ^2001bl) : lArnouts et arn2001bD : lYasuda et al.l 1I2OOII' ): ICapak et atl ||2004|) . Right panel 



Galaxy number counts in the i-band (not corrected for incompleteness) as compared to the FDF. The red dots show the average number 
counts of all COSMOS fields, whereas the green dots show the number counts as derived from the deep part of the FDF. The vertical 
dotted lines indicate the 50% completeness limits of COSMOS and FDF. The error bars show the Icr Poissonian errors. 




Figure 8. 1.5' x 1.5'Subaru i-band image. White circles mark the detected objects from the final catalogue, green circles mark objects 
eliminated from the final catalogues (saturated or corrupted magnitude, see text for details), whereas the blue circles mark false detections 
as detected on the negative image. It illustrates the fact, that there are many more detections in the negative image around bright objects 
(left panel) than in other areas of the image (right panel) . 
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Table 3. Galaxy number counts not corrected for incompleteness from COSMOS in the i (12 patches), Js (1 patch), H (15 patches), 
and K' (1 patch) bands, log and cri^g jv are given, where N is in units of mag~^deg~^. 
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1.721 


0.250 


1.492 


0.094 


2.022 


0.177 


16.5 


0.465 


0.307 


2.420 


0.112 


2.195 


0.042 


2.624 


0.088 


17.5 


1.465 


0.097 


2.721 


0.079 


2.738 


0.022 


2.925 


0.062 


18.5 


2.304 


0.036 


3.136 


0.049 


3.206 


0.013 


3.393 


0.036 


19.5 


3.127 


0.014 


3.525 


0.031 


3.631 


0.008 


3.790 


0.023 


20.5 


3.598 


0.008 


3.949 


0.019 


4.010 


0.005 


4.152 


0.015 


21.5 


3.996 


0.005 


4.277 


0.013 


4.236 


0.004 


4.283 


0.013 


22.5 


4.350 


0.003 


4.390 


0.011 


3.995 


0.005 


3.784 


0.023 


23.5 


4.677 


0.002 


3.705 


0.025 


2.486 


0.030 


1.244 


0.434 


24.5 


4.970 


0.001 


1.244 


0.434 


1.015 
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25.5 


5.103 


0.001 






0.170 
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0.001 
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120031 ) in constructing a multi-waveband catalogue in the 
COSMOS field. First we align by pixel all available filters 

(uCFHT, Bsubaru, Vsub aru; rSubarui ISubaru^ Zgubaru^ Kskpno, 

afterwards referred to as u, B, V, r, i, z, and K) to our H- 
band patches and derive i-selected catalogues. Please note 
that we had to compute a new astrometric solution for the 
public K band, as the solution given in the i mages was not 
sufiiciently accurate. The program SCAMP (jBertinl [iooel ') 
was applied for this procedure. Then we convolved all im- 
ages to the same seeing of 1.3"and ran SExtractor in the 
double image mode with the detection threshold of t = 2.5 
and n = 3 contiguous pixels (we detect on the original i- 
band image with a seeing of 0.95"). The 50 % completeness 
limits as well as the seeing of the different bands are listed 
in Table [l] The Table also compares the depth of the differ- 
ent bands with the depth in the PDF. Beside the u, z, and 
K bands, the COSMOS data set is roughly as deep as the 
PDF. In the u-band, z-band as well as in the K band the 
PDF is about 1 magnitude deeper. 

In the PDF we did the source detection with n = 3 con- 
tiguous pixels and a threshold of i = 1.7. This results in only 
a few false detections (less than 1 percent). For the COS- 
MOS data set we could not use the same detection threshold 
as this would result in too many false detections (measured 
on the negative image). We use a threshold of t = 2.5 in or- 
der to have false detections only at the 1 to 2 percent level. 
Depending on the patch, the false detections as measured 
on the negative image fluctuate between 0.6 percent and 2.6 
percent. Nevertheless this is not the real contamination rate, 
as a substantial fraction of these false objects detected on 
the negative image are distributed around bright objects. 
On the other hand the number density on the positive side 
does not increase around these bright objects. Therefore it 
is most likely, that these false detections are due to a reduc- 
tion problem triggering false detections only in the negative 
image. This can be seen in Fig. [S] where we show a region 
around relative bright objects together with the positive and 
negative detections. Taking this into account, we conclude 



that we have a contamination rate of about 1 percent or less 
in our catalogues. 

In total we detected about 300 000 objects in the 12 
patches (3 contiguous pixels and a detection threshold of 
2.5). These patches have a seeing of less then 1.3"(see above) 
in every band and reside in a region which has at least 90% 
of the maximum depth of each band. Moreover we exclude 
all objects (based on the weighting maps and SExtractor 
flags) with problems in the photometry (e.g. if some pixels 
of an object are saturated or the magnitude is corrupted) in 
at least one band in order to get a perfectly clean catalogue 
suitable to derive photometric redshifts. At this stage our 
clean i-selected catalogue comprises 293 377 objects. As we 
work on the individual patches and they slightly overlap, we 
analyse the objects and found that about 99% are unique 
entries. 

A comparison between the i-band number counts in the 
COSMOS field and in the literature as weU as with the PDF 
are shown in Pig. [T] There is a good agreement between the 
literature and COSMOS number counts up to the limiting 
magnitude. Moreover, the comparison with the PDF number 
counts shows an excellent agreement down to the faintest 
bins. Only at the very bright end the PDF number counts are 
slightly higher (most of the very bright objects are saturated 
in the COSMOS i-band and thus not taken into account), 
although not by more than Icr to 2a. The values of the galaxy 
number counts can be found in Table O 

In order to avoid contamination from stars, we rely on 
two sources of information: The star-galaxy classifier of the 
detection software SExtractor, and the goodness of fit of 
the photometric redshift code. We first exclude all bright 
{i < 22.5'") starlike objects (SExtractor star galaxy classi- 
fier > 0.97). Then we exclude all fainter objects whose best 
fitting stellar spectral energy distribution (SED) - according 
to the photometric redshift code - gives a better match to 
the fiux in the different passbands than any galaxy template 
(2 xltar < xlaiaxy)- These objects are subsequently flagged 
as stars and removed from our galaxy catalogue. In order 
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to test the accuracy of our procedure we further inspect by 
eye (in one patch) on the public ACS data if the objects 
flagged as stars are extended. It turned out that about half 
of these objects are extended in the ACS dat a and only half 
are po int-like. Although this approach (see iGabasch et al.l 
l2004al ) of excluding stars works very well in the FDF where 
the seeing of the detection image is only 0.55", it is much 
less effective in the COSMOS i band with a seeing of 0.95". 

In total 4803 (1.6 %) objects were classified as stars and 
removed from our sample (see also Fig. [TS] for the photo- 
metric redshift distribution of all starlike objects excluded 
from the galaxy catalogue). Please note that most of the 
bright point-like objects are saturated in one of the bands 
and already removed in the first cleaning step (see above). 
Therefore our final i-band selected galaxy catalogue com- 
prises 288 574 objects. 



5 SPECTROSCOPIC REDSHIFTS 

Photometric redshifts need to be calibrated by spectroscopic 
redshifts for objects covering a wide range of galaxy types 
and redshifts. H ence, we selected sp ectroscopic data of the 
COSMOS field (|ScoviUe et al.ll2007l ) f rom the ESO archive 
taken with the VIMOS spectrograph l|Le Fevre et aLlbOOal ) 
at the ESQ VLT in the context of the zCOSMOS project 
l|Lillv et al.ll2007l ). zCOSMOS is a Large ESO Programme 
with 600 hr of observing time aiming at the characteri- 
sation of the distribution and properties of galaxies out 
to redshifts of z ~ 3. The project is divided into two 
parts. zCOSMOS-bright focuses on relatively bright galax- 
ies with 7ab < 22.5 at redshifts 0.1 < z < 1.2. In or- 
der to include such galaxies in our spectroscopic control 
sample we selected the VIMOS mask 'zCOSMOS_4-4' with 
312 slits centred at 10h00m28.20s -|-02°15'45.0". This mask 
was observed 1.5 h with the red low-resolution LR_red grism 
{R ~ 210), covering a wavelength range from 5500 to 9500 A. 
zC OSMOS-deep aims at galaxies with 1.5 < z < 2.5. For 
this lLiUv et al. (2007) selected targets using the BzK crite- 
ria of iDaddi et al. (200 4') arid the UGR 'BX' and 'BM' se- 
lection of lSteidel et al ] |2004l l. We obtained a representative 
sample of such galaxies, reducing the VIMOS mask 'zCOS- 
MOS_55_faint' with 239 slits at 10h00m27.67s -h02°10'23.0", 
observed 4.5 h with the blue low-resolution LR_blue grism 
[R ~ 180), covering a wavelength range from 3700 to 6700 A. 

The sp ectra of both VIMOS masks were reduced us- 
ing VIPGI IScodeggio et al.ll2005l ). This software package 
is designed to reduce the multi-object VIMOS spectra in 
a quite automatic way. The reduction was performed us- 
ing essentially standard methods. The jittered sequences of 
individual exposures (five different positions) allowed to effi- 
ciently correct for fringing and to obtain a high-quality sky 
subtraction. The one-dimensional compo site spectra w ere 
extracted by means of the S/N-optimised iHornd (|l986l ) al- 
gorithm. The final spectra are fiux calibrated and corrected 
for atmospheric absorption bands. 

Redshifts were derive d using the cro ss-correlation based 
algorithm described bv iNoU et al.l (|2004h . In order to obtain 
the redshift and a rough spectra l type, we used a sequence 
of six empirical templates (see INoII et al]l2004l . for more 
details) essentially differing in their UV-to-optical flux ra- 
tio, but also showing different strengths of nebular emission 
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Figure 9. Redshift distribution of the COSMOS spectroscopic 
sample for (a) all extragalactic objects with trustworthy redshifts, 
(b) the corresponding objects of the bright-galaxies mask, (c) 
the corresponding objects of the zCOSMOS-deep mask, (d) - (h) 
galaxies of different object type, and (i) the quasars. The types 
range from I = ellipticals ( or passively evolving) to V = extreme 
star bursts (see INoU et aLll2OO40 . In diagram (h) galaxies of type 
V (white bars) as well as type VI (strong Lya emitters, black 
bars) are plotted. The redshift resolution is Az = 0.05. 



lines. For the comparison to photometric redshifts we con- 
sider galaxies with trustworthy (> 90% confidence) spectro- 
scopic redshifts only. Due to low S/N (particular at high 
redshift) and the limited wavelength ranges of the spec- 
tra certain redshifts could only be derived for about half 
of the objects, i.e. 272 of 551 spectra. For the zCOSMOS- 
deep mask the success rate (33%) is significantly lower than 
for the bright-galaxy mask (62%). The total sample of 272 
objects with trustworthy redshifts includes 49 stars and 12 
quasars. Hence, the final sample of galaxies comprises 211 
objects. The majority of them (147 or 70%) belongs to the 
bright-galaxy sample. On the other hand, the subsample of 
36 galaxies with redshifts 1.5 < z < 3 only consists of ob- 
jects from the zCOSMOS-deep sample. These high-redshift 
galaxies (quasars excluded) correspond to 45% of the ob- 
jects identified in the deep sample. Fig. [5] shows the red- 
shift distribution of the full spectroscopic control sample, 
the 'zCOSMOS_4-4' mask, the 'zCOSMOS_55_faint' mask, 
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the different galaxy types defined by iNoll et all (|2004 I = 
ellipticals, V = extreme starbursts, VI = strong Lya emit- 
ters), and the quasars. 

After cross-correlating (by visual inspection) the 211 
spectroscopic objects with the final i-band selected cata- 
logue in the 12 patches we end up with a final sample of 
162 spectroscopic redshifts used to calibrate the photometric 
redshifts. As the spectroscopic redshift sample was merely 
used to calibrate the photometric redshift and to measure 
the accuracy of the photometric redshift code, we decided to 
use only objects with an unique cross identification. There- 
fore we excluded all objects where the cross-matching be- 
tween the photometric and spectroscopic catalogue was not 
unique, where the objects of the photometric catalogue show 
saturation in one (or more) image and where the spectro- 
scopic distance determination was not very reliable (e.g. low 
S/N ratio). 



6 PHOTOMETRIC REDSHIFTS 

A summary of the photometric redshift technique used 
to derive the dista nces to the galaxie s can b e found in 
iBender et all (120011 ) and lGabasch et al.l l|2004al '). Before de- 
riving the photometric redshifts we checked and fine-tuned 
the calibration of our photometric zeropoints by means of 
colour-colour plots of stars. We compared the colours of 
st ars with the co lours of stellar templates from the library 
of lPicklesI l|l998h converted to the COSMOS filter system. 
In general, corrections to the photometric zeropoints of only 
a few hundredth of a magnitude were needed to obtain an 
good match to the stars and best results for the photometric 
redshifts (if compared to the spectroscopic ones). Only in the 
u-band and in the K-band, the correction were in the order 
of a few tenths of a magnitudes. A comparison between the 
reduced KPNO K-band and our K-band (both convolved to 
the same seeing of 1.3") in patch 01c showed, that although 
the total magnitudes agreed very well, the fixed aperture 
magnitudes (especially of relatively faint sources) differ sys- 
tematically by a few tenths of a magnitude. Therefore we de- 
cided to correct the KPNO zeropoint by matching the faint 
sources in the KPNO image with those in our own field (in 
the fixed aperture used to derive the photometric redshifts) . 
Moreover when calculating the photometric redshifts, we ar- 
tificially increased the magnitude errors in the K-band by 
0.25™ (added in quadrature to the SExtractor errors) to re- 
duce the relative weight of this slightly problematic band. 
Therefore, we rely mostly on the accurate photometry of the 
NIR H band. 

In order to avoid contamination from close-by objects, 
we derived object fiuxes for a fixed aperture of 2.0" (1.5x 
seeing) from images which had been convolved to the same 
point spread function (PSF; 1.3"). A redshift probability 
function P(z) was then determined for each object by match- 
ing the object's fiuxes to a set of 29 template spectra red- 
shifted between z = G and z = 10 and covering a wide range 
of ages and star-formation histories. 

In Fig. [To] (left panel) we compare 162 high quality 
galaxy spectroscopic redshifts with the photometric red- 
shifts. Although there is a good agreement in the redshift 
range between z ~ 0.2 and z ~ 1.2, it is clear from Fig. 1101 
that there is a degeneracy between high redshift (2 ~ 2.5) 



and low redshift [z ~ 0.2) objects (10 catastrophic outliers 
with /S.z/{zspEc + 1) > 0.2). This degeneracy stems from the 
relatively red u-band. In Fig.[T2]we show the redshift proba- 
bility function as well as the SED fits to the observed fiux of 
the spectroscopic object 000932543. Although the spectro- 
scopic redshift is Zspec = 0.093, the best fitting photometric 
redshift is Zphot = 2.72. On the other hand, there is also a 
low redshift peak around Zphot = 0.18 in the redshift prob- 
ability function (but with a lower probability). Moreover 
Fig. m] also shows that both, the high redshift as well as the 
low redshift solutions are hard to disentangle as long as no 
information in the UV is available (as they differ mainly in 
the UV). Therefore we decided to include in the determina- 
tion of the photometric redshifts also the GALEX FUV and 
NUV bandfl 

As we do not want to convolve all the images to a seeing 
of 5" (GALEX PSF), we decided to use another approach to 
include the UV fluxes in our photometric redshift estima- 
tion. Similar to the optical and NIR bands we used a fixed 
aperture of ~ 1.5x PSF, i.e. 7.5". As there were no obvious 
features in the colour-colour plots of stars including the UV 
bands, we could not fine-tune the calibration of the zero- 
points by means of colour-colour plots of stars. Therefore 
we optimised the zeropoints by using the SED fits of our 
galaxies with very good photometric redshifts (if compared 
with the spectroscopic ones) . Please note that this approach 
can not derive accurate UV fluxes, but gives only a very 
rough flux estimation in the two UV bands. Nevertheless, it 
is now possible to break the degeneracy between the high 
redshift and low redshift solution. This can be best seen in 
Fig. [T2]where we show one of the catastrophic outliers. Only 
by including the NUV and FUV fiuxes (right panel) we are 
able to obtain a photometric redshift oi z = 0.08, hence very 
close to the spectroscopic redshift of Zspec ~ 0.093. This ap- 
proach drastically reduces the number of our catastrophic 
photometric redshift outliers (see Fig. I10|l . 

In Fig. [TO] (right panel) we compare the final photo- 
metric and spectroscopic redshifts of the 162 galaxies. The 
agreement is very good and we have only 3 catastrophic out- 
liers. The right panel of Fig. llll shows the distribution of the 
redshift errors. It is nearly Gaussian and scatters around 
zero with an rms error of Az/{zspsc + 1) ~ 0.035. Fig. [TT] 
(left panel) presents the distribution of the best fitting 
templates and photometric redshifts for all the objects. The 
median value of the reduced is 1.5 and demonstrates that 
the galaxy templates describe the vast majority of galaxies 
very well. Please note that although the photometric redshift 
accuracy of a single object does not considerably improve 
by adding our H-band to the publically available data, the 
number of catastrophic outliers decreases by nearly a factor 
of 2. Moreover, the H-band was also very useful to find and 
address the problems in the public K-band as mentioned be- 
fore. 

The galaxy redshift histogram of all objects in the different 
patches is shown in Fig. 1131 The mean galaxy distribution 
can be very well des cribed by Equation (O introduced by 
iBrainerd" et all (|l996l ). 



* The data were taken from: 
|http://irsa.lpac! caltech.edu/data/COSMOS/ 1 
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Figure 10. Comparison of spectroscopic and photometric redshifts in the COSMOS field (162 galaxies). Left panel: u band to K band 
are used to derive the photometric redshifts. Right panel: u band to K band as well as the GALEX FUV and NUV bands are used to 
derive the photometric redshifts. 
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Figure 11. Left panel: Histogram of the reduced for all galaxies in the COSMOS field as obtained for the best fitting template and 
redshift. The dotted vertical line indicates the median reduced of 1-5- Right panel: Histogram of the photometric redshift errors. The 
error distribution can be approximated by a Gaussian centred at 0.007 with an rms of 0.035 (solid line). 



Pz{z) — Const X 



exp 



(-(^Ao) 



(7) 



where Const, zo =< z > and (3 are free parameters 

with < z > being the first moment of the distribution and F 
the Gamma function. The best fitting values are: Const — 
6206, zo = 0.107, and f3 = 0.611. 



Please note that if we analyse the galaxy photometric 
redshift histogram with a binning of Az = 0.1 there are 
three clearly visible peaks below redshift of z — 2: one at 
Zphot = [0.6, 0.7], one at Zphot = [0.9, 1.0], and one at Zphot = 
[1.7, 1.8]. Interestingly, we also find peaks in the spectro- 
scopic redshift histogram (see also Fig. [UJ with at least 10 
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Figure 12. Redshift probability distribution (upper part) and photometric redshift fit (lower part) of the spectroscopic object 000932543 
with Zspec = 0.093. The difi'erent SEDs are colour coded. The open circles (lower part) represent the SED integrated within the various 
filter transmission curves, whereas the black dots represent the measured fluxes. The redshift for the best fitting SED (^i) as well as for 
the second best fit (22) are also given. Left panel: Only 8 bands (u-band to K-band) are used to derive the photometric redshift yielding 
a wrong Zphot of 2.72. Right panel: 10 bands (FUV, NUV, u-band to K-band) arc used to derive the photometric redshift yielding a 
Zpfiot of 0.08 very close to the spectroscopic redshift. 
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Figure 13. Left panel: Redshift number distribution of all galaxies in the 12 patches. The mean number d i stribu tion is shown as blue 
dots. The blue solid line represents a fit to the mean galaxy distribution using the approach of lBrainerd et al] lll996l) (see text for details). 
Right panel: Redshift distribution of objects classified as stars (see text for details) and eliminated from the final galaxy catalogues. 
Please note the different scaling by a factor of 20 in the left and right panel. 
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Figure 14. Absolute UV-band magnitudes of galaxies in the 
COSMOS 06a patch (black dots) as a function of redshift. The 
green dots are derived in the FDF whereas the red dots represent 
COSMOS objects classified as stars if a more conservative crite- 
rion in separating stars from galaxies (Xstar < X^alaxy ^ 
see text for details) is applied. 
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7 UV LUMINOSITY FUNCTION AND 
DENSITY 

In this section we show the UV luminosity function (LP) at 
1500 A as derived from our deep i- selected catalogue an d 
compare it to the LPs in the FDP (jCabasch et al.l l2004al ) . 
Note that in this paper we are not aiming to give a complete 
analysis of the LP evolution. A detailed analysis of the LPs 
derived in different passbands as well as the star-formation 
rate (SFR) together with an elaborate comparison to the 
recent literature will be presented in future papers. Here we 
limit ourselves to a one-to-one comparison with the PDP 
luminosity functions (mostly as a consistency check to our 
earlier work). 

To derive the absolute UV band magnitude we use the 
best fitting SED as determined by the photometric red- 
shift code. Since the photometric redshift code works with 
aperture fluxes, we only need to correct to total luminosi- 
ties by applying an object dependent scale factor. Por this 
scale factor we used the ratio of the I-band aperture flux to 
the total flux as provided by SExtractor (MAG_APER and 
MAG_AUTO). As the SED fits all observed- frame passbands 
simultaneously, possible systematic errors which could be in- 
troduced by using K-correcti ons applied to a single observed 
magnitude are reduced (see iGabasch et al.ll2004al . for more 
details) . 

As an example, we plot in Pig.[T4]the absolute UV-band 



magnitudes against the photometric redshifts of the objects 
in the COSMOS patch 06a. Moreover we also show the abso- 
lute UV-band magnitudes as derived in the PDP. Both fields 
agree very well in their magnitude distribution, although 
there are a few relatively bright objects in the COSMOS 
patch (about 5 times the area of the PDP) not seen in the 
PDP distribution. To check if those objects could be stars 
misclassified as galaxies by our star-galaxy separation crite- 
rion, we decided to use a more conservative criterion for sep- 
arating stars from galaxies. We changed our criterion from 

2 Xitar < Xlalaxy tO X^star < X^alaxy (see above) . As CaU be 

seen in Pig. [T?] even this conservative criterion does not re- 
move a substantial fraction of these relatively bright objects. 
Note that because of the larger seeing (0.95") compared to 
the PDP (0.55"), there may be more blended objects in the 
catalogue. 

The luminosity function is computed by dividing the 
number of galaxies in each magnitude bin by the volume 
Vbin of the redshift interval. To account for the fact that 
some fainter galaxies are not visi ble in the whol e survey vol- 
ume, we performed a V/Vmax (|Schmidtl Il96 8^) correction. 
The errors of the LPs were calculated by means of Monte- 
Carlo simulations and include the photometric redshift error 
of every single galaxy, as well as the statistical error (Pois- 
sonian error). To derive precise Schechter parameters, we 
limited our analysis of the LP to the magnitude bin where 
the V/Vmax correction is negligible (red dots in Pig. I15p . 
We also show the uncorrected LP in the various plots as 
open circles. We did not assume any evolution of the galax- 
ies within the single redshift bins. The redshift intervals are 
approximately the same size in ln(l -I- z), and most of the 
results we are going to discuss are based on 1000 - 4000 
galaxies per redshift bin and per patch. 

In Pig. [T^ we present the UV luminosity functions at 
1500 A (we evaluate the luminosity function in the rect- 
angular filter at 1500 ± 100 A). The filled (open) symbols 
denote the luminosity function with (without) completeness 
correction in the different patches. We also show the V/Vmax 
corrected mean LPs in the CO SMOS field as well as the 
PDP LPs l|Gabasch et al.ll2004al '). The solid red hues show 
the Schechter function fitted to the luminosity function (we 
used a fixed slope of a = —1.07 as found in the PDP). The 
best fitting Schechter parameter, the redshift binning as well 
as the reduced a-re also listed. 

It is obvious from the figure that there is strong evo- 
lution in characteristic luminosity and number density be- 
tween redshifts 0.6 and 4.5. This can be best seen in Pig. 1161 
where we show the redshift evolution of M* and (p* as derived 
from the Schechter functions fitted to the LPs of Pig. 1151 
Please note that we used only magnitude bins with a V/Vmax 
correction of about unity. Moreover we exclude also very 
bright magnitude bins as there might be a contamination 
by spurious detections, stars, AGNs or blended objects. The 
magnitude bins used for deriving the best fitting Schechter 
parameters are shown in Pig. [15] as red dots. Nevertheless, 
including also the very bright bins changes the best fitting 
values only negligibly. We find a substantial brightening of 
M* and a decrease of </>* with redshift: from (z) ~ 0.5 to 
(z) ~ 4.5 the characteristic magnitude increases by about 3 
magnitudes, whereas the characteristic density decreases by 
about 80 - 90%. Note that our result s do not disa gree with 
recent findings of e.g. iBouwens et al.l (|2006l . |2007| ) who are 
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Figure 15. Luminosity functions at 1500 A from low redshift ((2) = 0.3, upper left panel) to high redshift ((2) = 5.5, lower right panel). 
The filled (open) black symbols show the luminosity function corrected (uncorrected) for V/Vmax in the various patches. The red and the 
blue dots represent the mean LF in the field. The green squares represent the LFs as derived in the FDF. The fitted Schechter functions 
for a fixed slope a are shown as solid red lines (only the red dots are used to derive the best fitting Schechter functions). Note that we 
only fit the luminosity functions to (2) = 4.5. The parameters of the Schechter functions can be found in Table [5] The Schechter fit for 
redshift (2) = 0.6 is indicated as a dashed black line in all panels. 



tracing the UV LF from z; ~ 3 to 2 ~ 6 and find a fainten- 
ing of M* (with increasing redshift) and nearly no density 
evolution in this redshift range since we limit our ana l ysis to 
redshifts of 2: < 4.5. As discussed in iBouwens et al.l (|2007l ) 
it is plausible that there is a turnover in the LF evolution 
at redshift of z ~ 4. 

The best fitting Schechter parameters and their la errors in 
the COSMOS field are summarised in Tabled In Fig.[T6lwe 
also show the redshift evolution of M* and (p* as derived in 



the FDF bv I Gabasch et alj ()2004al ). Although in the COS- 
MOS field the characteristic magnitude is about 0.25 mag 
brighter in most of the redshift bins if compared to the FDF 
results, the characteristic density is slightly lower. This re- 
sults in a very similar UV luminosity density (LD; the inte- 
grated fight emitted by the galaxies) in the COSMOS field 
and in the FDF, as can be seen in Fig. 1171 
The LD at a given redshift is derived by summing the 
completeness-corrected (using a V/Vmax correction) lumi- 
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Figure 16. Redshift evolution of M* (left panel) and 0* (right panel) for the UV band at 1500 A. The red dots are derived from the 
mean luminosity function in the COSMOS field (Schechter fit to the red dots in Fig. I15I I whereas the green dots stem from the FDF. 



Table 5. Schechter function fit at 1500 A 
redshift interval 



M* (mag) 



(Mpc-3) 



a (fixed) 



0.41 - 0.71 -18.20 +0.09 -0.09 1.12e-02 +7.61e-04 -7.61e-04 -1.07 

0.71 - 1.11 -19.37 +0.05 -0.05 5.66e-03 +1.97e-04 -1.97e-04 -1.07 

1.11 - 1.71 -19.66 +0.04 -0.04 5.38e-03 +1.97e-04 -1.69e-04 -1.07 

1.71 - 2.15 -20.28 +0.08 -0.08 3.77e-03 +3.94e-04 -3.38e-04 -1.07 

2.15 - 2.91 -20.70 +0.09 -0.09 2.90e-03 +3.94e-04 -3.38e-04 -1.07 

2.91 - 4.01 -21.01 +0.08 -0.09 1.58e-03 +1.41e-04 -1.41e-04 -1.07 

4.01 - 5.01 -21.32 +0.14 -0.13 5.35e-04 +5.63e-05 -5.63o-05 -1.07 



nosity of every single galaxy up to the abso lute magnitude 
limit. C ontr ary to the procedure d escribed in lGabasch et al] 
l|2004bl ) and lGabasch et al.l (|2006l ). we do not apply a further 
correction (to zero galaxy luminosity), to take the missing 
contribution to the LD of the fainter galaxies into account. 
This further correction should be done only after a very care- 
ful analysis of the LFs (as it requires an extrapolation of the 
LF at the faint end) and is postponed to a future paper. As 
the limiting i-band magnitudes of the COSMOS field and 
the FDF are very similar, we decided to integrate the LF 
in both fields down to the same absolute magnitude instead 
(see Table [6]). This approach allows us to directly compare 
the LD of both fields without relying on any extrapolation. 

In Fig. [l7]we show the 1500A luminosity densities of the 
single COSMOS patches, the mean COSMOS LD as weU as 
the LD derived in the FDF. The FDF and COSMOS lu- 
minosity densities are integrated down to the faint-end lim- 
iting magnitudes given in Table |6l Moreover we also show 
the COSMOS luminosity densities integrated between the 
faint-end and a bright-end limiting magnitude (see Table [B]). 
The bright-end cut excludes in total 80 (0.03%) bright ob- 
jects in the COSMOS field. It was derived by comparing the 
Schechter-fit to the LF (the excluded LF bins are shown as 



blue dots at the bright end in Fig llSp as well as by com- 
paring the absolute UV-band magnitudes as a function of 
redshift with those in the FDF (see Fig. [14] for one patch). 
Fig. [T7] nicely shows that the UV LD and its redshift evolu- 
tion in the FDF (based on a very small field if compared to 
the COSMOS field) agrees quite well with the result in the 
COSMOS field. Even though for redshifts below z ~ 2.5 the 
UV LD in the FDF is systematically lower than the value in 
the COSMOS field the deviation is in the order of la. The 
values of the mean COSMOS LD together with the error 
bars are listed in Table [S] 



8 SUMMARY AND CONCLUSIONS 

In this paper we present the data acquisition and reduction 
of NIR Js, H, and K' bands in the COSMOS field. We de- 
scribe a 2-pass reduction pipeline to reduce NIR data. The 
2-pass pipeline is optimised to avoid flat-field errors intro- 
duced if the latter are constructed from science exposures. 
Moreover we present and implement a method to stack im- 
ages of different quality resulting in an optimal S/N ratio 
for faint sky dominated point sources. The Js and K' band 
cover an area of about 200D' (1 patch) whereas the H band 
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Figure 17. The 1500A luminosity densities of the single COSMOS patches (black dots), the mean COSMOS LD (red dots) as well as 
the LD derived in the FDF (green dots) . Left panel; The luminosity densities are integrated between the faint-end limiting magnitudes 
given in Table[6]but no limiting magnitude cut was used for the bright-end. Right panel: The luminosity densities are integrated between 
the faint-end and the bright-end limiting magnitudes given in Table [6] The values of the mean COSMOS LD together with the error 
bars are listed in Table [6] 



Table 6. The mean COSMOS LD at 1500 A for the different red- 
shift bins. The luminosity densities are derived within the faint- 
end and the bright-end limiting magnitudes given in the last col- 
umn. 



redshift 



luminosity density 



magnitude 







(W Hz-i Mpc-3) 


ran 


?e 


0.15 


- 0.41 


7.557e+18 ± 7.243e-|-17 


-14. - 


-22. 


0.41 


- 0.71 


9.640e+18 ± 7.219e-|-17 


-14. - 


-22. 


0.71 


- 1.11 


1.536e+19 ± 7.329e-|-17 


-15. - 


-23. 


1.11 


- 1.71 


1.750e+19 ± 6.290e-|-17 


-16. - 


-23. 


1.71 


- 2.15 


2.097e+19 ± 1.719o+18 


-17. - 


-23. 


2.15 


- 2.91 


2.340e+19 ± 2.116c-M8 


-18. - 


-24. 


2.91 


- 4.01 


1.485e-fl9 ± 6.931e-|-17 


-19. - 


-24. 


4.01 


- 5.01 


7.809e-f 18 ± 4.728e-|-17 


-19. - 


-24. 


5.01 


- 6.01 


6.052e-f 18 ± 8.592e-(-17 


-19. - 


-25. 



covers about 0.85D° (15 patches) in total. The 50% com- 
pleteness limits are 22.67, ~ 21.9, and 21.76 in the Js, H, 
and K' band, respectively. The number counts of all NIR 
bands nicely agree with the number counts taken from lit- 
erature. 

Furthermore we present a deep and homogeneous i-band 
selected multi- waveband catalogue in the COSMOS field by 
combining publicly available u, B, V, r, i, z, and K bands 
with the H band. The clean catalogue with a formal 50% 
completeness limit for point sources of i ~ 26.7 comprises 
about 290 000 galaxies with information in 8 passbands and 
covers an area of about 0.7D° (12 patches). We exclude all 
objects with corrupted magnitudes in only one of the filters 
from the catalogue in order to have a catalogue as homoge- 
neous as possible. 



Photometric redshifts for all objects are derived and a 
comparison with 162 spectroscopic redshifts in the redshift 
range < 2: < 3 shows that the achieved accuracy of the 
photometric redshifts is Az/(zspec + 1) < 0.035 with only 
~ 2% outliers. Please note that in order to break the de- 
generacy between high redshift and low redshift solutions 
we included also the GALEX FUV and NUV fihers in the 
photometric redshift estimation which considerably reduced 
the number of outliers. 

The multi-waveband catalogue including the 
photometric redshift information is made pub- 
licly available. The data can be downloaded from 
\protect\vrulewidthOp t|http://www.mpe.mpg.de^gabasch/COSMOS/| 

We derive absolute UV magnitudes and a comparison 
in a magnitude-redshift diagram with the FDF shows good 
agreement. Moreover we investigate the evolution of the lu- 
minosity function evaluated in the restframe UV (1500 A). 
We find a substantial brightening of M* and a decrease of 
(j)* with redshift: from {z) ~ 0.5 to {z) ~ 4.5 the character- 
istic magnitude increases by about 3 magnitudes, whereas 
the characteristic density decreases by about 80 - 90%. 

We compare the redshift evolution of the UV luminosity 
density in the COSMOS field and the FDF up to a redshift 
of 2: ~ 5. Below a redshift of z ~ 2.5 the mean UV luminosity 
density in COSMOS is systematically higher by about Icr if 
compared to the FDF. At 2.5 < 2 < 5 both UV luminosity 
densities agree very well. 

It is worth noting the remarkably good agreement be- 
tween the UV LF as well as the UV LD despite the fact, 
that the FDF is about 60 times smaller than the COSMOS 
field analysed here. 
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